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Abstract 


The Environmentally Responsible Aviation project seeks to accomplish the 
simultaneous reduction of fuel burn, noise, and emissions. A project at NASA 
Dryden Flight Research Center is contributing to ERAs goals by exploring the 
practical application of real-time trim configuration optimization for 
enhanced performance and reduced fuel consumption. This peak-seeking 
control approach is based on Newton-Raphson algorithm using a time- 
varying Kalman filter to estimate the gradient of the performance function. In 
real-time operation, deflection of symmetric ailerons, trailing-edge flaps, and 
leading-edge flaps of a modified F-18 are directly optimized, and the 
horizontal stabilators and angle of attack are indirectly optimized. 
Preliminary results from three research flights are presented herein. The 
optimization system found a trim configuration that required approximately 
3.5% less fuel flow than the baseline trim at the given flight condition. The 
algorithm consistently rediscovered the solution from several initial 
conditions. These preliminary results show the algorithm has good 
performance and is expected to show similar results at other flight conditions 
and aircraft configurations. 


Agenda 

• Background 

— Motivation 

- Previous Approaches 

- New Approach 

• Precursor Open-Loop Experiment 
Flight Results (1 Flight) 


• Closed-Loop Algorithm Preliminary 
Flight Results (2 Flights) 
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NASA Goals 


• This work directly supports the Environmentally Responsible 
Aircraft (ERA) Project’s goal of reducing fuel burn. 


• Aeronautics Research Mission Directorate 

- Integrated System Research Program 

• Environmentally Responsible Aviation 

- Airframe Technologies 

>> 2.2 Flight Dynamics and Controls 

• 2.2.4 Intelligent Controls 


Motivation 


• Multiple longitudinal effectors 
for trim 

- Traditionally horizontal tail 
incidence angle or elevator. 

- But also: Symmetric ailerons, 
flaps, leading-edge devices, thrust 
vectoring, pump fuel fore/aft for 
c.g. control, etc. 

• Is there an alternative, lower- 
drag trim solution? 

• Can we adjust to variations 
between: 

- Aircraft? 

- Configurations? 

- Flight conditions? 


Section Lift Coefficient 


Initial Trim 




Modified Trim 


Previous Research Results 



Adaptive Performance Optimization 

Patent 5,908,176 

Gilyard's L-1011 flight test results in 1999: 

“Optimizing the symmetric outboard aileron 
position realizes a drag reduction of 2-3 drag 
counts (approximately 1 percent)." 
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Symmetric aileron, deg 

Figure 5. Variation of incremental drag with symmetric 
outboard aileron deflection for a two-sided raised-cosine 
maneuver (data filtered, period = 150 sec, 

a = 0.000125). 


Flight Test of an Adaptive Configuration Optimization 
System for Transport Aircraft 

Gilyard, Glenn B.; Georgie, Jennifer; Barnicki, Joseph S. 
Dryden Flight Research Center, 1999. 



http://hdl.handle.net/2060/1999001943S 




Dryden FI i ght Research Center EC97-44077-3 Photographed May 1 997 
This modified Lockheed L-1011 TriStar operated by Orbital 
Sciences Corp. lifts off on its first flight in NASA’s Adaptive 
Performance Optimization study. (NASA/Tony Landis) 


Performance Improvement Package (PIP) for 111 


Figure 2: Drooped aileron 

Boeing engineers determined that a 2-degree aileron droop was optimal for flight performance. 



Boeing, United Teaming To Improve Fuel Efficiency. 

The International Business Times (3/23, Francheska) reports, "Boeing 
and United Continental Holdings, Inc. has entered into an agreement to 
modify United Airlines' 777 fleet with a Performance Improvement 
Package with the aim of achieving greater fuel efficiency and reduced 
emissions." The upgrade "improves the airplane's aerodynamics 
through a software change to enable a drooped aileron, a ram air 
system improvement and the installation of improved wing vortex 
generators." If gas costs $100 per barrel, the program is expected to 
save each plane $200,000 a year in gas costs. 

http://au.ibtimes.com/articles/125677/20110323/boeing-united-airlines-fleet-airline-airplane-fuel-777.htrn 



Delivering Fuel and Emissions Savings for the 777 
By Ken Thomson, and E. Terry Schulze 

http://www.boeing.com/commercial/aeromagazine/articles/qtr 03 09/pdfs/AERQ Q309 article02.pdf 


Boeing Trailing Edge Variable Camber (TEVC) System 



TEVC System on 787: 

"The TEVC cleverly articulates the trailing edge of the 
flaps in various cruise conditions to help reduce drag." 

Guy Norris, Aviation Week in 2010 

http://www.aviationweek.com/aw/blogs/commerdal aviation/ThingsWithWings/index.jsp?plckController=Blog&plckScri 
pt=blogScript&plckElementId=blogDest&plckBlogPage=BlogViewPost&plckPostId=Blog%3A7a78f54e-b3dd-4fa6- 

ae6e-dff2ffd7bdbbPost%3A57b52637-d9a6-4589-bl74-7aece76b0c46 


However, the new 747-8 does not have it: 

"Boeing also acknowledged that the Trailing Edge 
Variable Camber (TEVC), one of the unique features 
on the 787 Dreamliner, does not feature on the 747-8 
family. 

'The 747-8 does not have TEVC. The 747-8 Program did 
study adding TEVC but found the performance benefit 
was not as great as on the 787. The 747-8 does not 
have active gust suppression control laws. The 747-8 
has better response to wind gust due to its physical 
size. However, the 747-8 will have active control laws 
designed to improve the ride quality of the airplane.'" 

Daniel Tsang, Aspire Aviation in 2010. 

http://www.aspireaviation.com/2010/06/21/747-8f-flight-test-progressing/ 




Case for In-Flight Trim Optimization 




Aircraft design cycles are long. 
Fleets are aging. 


Real-time optimization can adapt to 
conditions unforeseeable during the 
design phase. 


New Approach 

Goal: 

Reduce fuel burn in cruise flight. 

Technical Challenge: 

- Implement an algorithm to change trim allocation in flight for minimum fuel 
consumption. 

- Evaluate if the algorithm is well-behaved in realistic conditions. 

General Approach: 

- Move effectors to minimize drag to reduce fuel flow. 

- Maintain trim cruise speed and altitude. 

- Peak-seeking control, steepest descent algorithm. 

- Estimate performance function gradient with time-varying Kalman filter. 

- Simultaneously optimize the position of multiple effectors. 

- Kalman filter addresses noise issues directly. 

- New method not reliant on models. 

- Controller seeks the minimum fuel flow solution. 


Peak-Seeking Control 

• Given: f 15 

I 10 

- A performance measurement, fuel flow, * 5 

that is a function of surface positions | o 

• The minimum-cost (blue) combination of 5 

surface positions (x,y,z) is unknown 

• This is called the Performance Function 

- Measurements of surface positions and fuel 
flow are noisy. 

• Find: 

- Minimum of the performance function, in flight 

• Assumptions: 

- Performance function has a single minimum 

- Measureable surface positions and fuel flow 

- Gaussian distributed noise 

- Plant is stable and controllable (inner loop 
control design treated as separate problem) 


Notional Performance Function 
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Symmetric Aileron (deg) 

Fuel Flow (color) due to 
Surface Deflections (x,y) 



Fuel Flow (color) due to 
Surface Deflections (x,y,z) 




Performance Measurements 


Simplified Introduction to Approach: Single Effector 



Effector Position, x 

(Commanded by Peak-Seeking Controller) 


Performance Measurements 


Simplified Introduction to Approach: Single Effector 
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Simplified Introduction to Approach: Single Effector 
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Performance Measurements 


Simplified Introduction to Approach: Single Effector 



Effector Position, x 

(Commanded by Peak-Seeking Controller) 


Technical Approach: Architecture 


Peak seeking control approach based on work by Ryan and Speyer: 

Ryan, J.J. and Speyer, J.L., "Peak-Seeking Control Using Gradient and Hessian Estimates" 
Proceedings of the 2010 American Control Conference, June 30-July 2, 2010, pp. 611-616. 

http://hdl.handle.net/2060/20100024511 



Taylor series expansion of the performance function: 

f{x) * f(x k ) + bl(x-x k ) + O jx. - x. k ) 

Gradient Higher order terms 

Estimate 
(transpose) 

19 


Technical Approach: Performance Function 


Change in PF: 


Change in surface positions: 


AA = /Ok-i) - / Ok) 



y y 

k ~ A /c— 1 


Assuming PF can be treated as linear at any control surface position: 

A/k = blAx k 


In this 2D example, we have two 
groups of effectors, symmetric 
ailerons and TEF: 


Ax 2t 



k—1 


k - 1 


- X 


-X 2 


1 k 

k 


Rewritten: 

a/i = 





(Ryan 2010) 


This has been extended to 3D, and can be extended to more effectors. 
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Technical Approach: Kalman Filter States 



Since the gradient may change with surface positions and measurements will be noisy, 
a Kalman filter is an appropriate choice of estimator. 


Kalman filter states chosen to be: 


6c = 


b 

b- 


k 


2 k 


The measurement equation of the linear 
time-varying Kalman filter takes the form: 


A//C — H k Zk + v k 



process and the linear time-varying Kalman filter process 
equation is then given by: 


6c+ 1 - + w k 


(Ryan 2010) 
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Linear Time- Varying Kalman Filter 


Current Kalman gain: 

Current KF state estimate (gradient of PF): 
Current state covariance matrix: 

Predicted KF state estimate (gradient of PF): 
Predicted state covariance matrix: 


K = P k Hl{H k P k Hl + V k ) 

<* = 4 + K(Af k - h * 4 ) 

P k = 0 - KH k )P k 

Cfc+l Zk 

?k+l = ?k +W k 


The noise variances Vk and Wk are used as 
tuning parameters influencing the filters 
performance. 


(Ryan 2010) 
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Flight Research Implementation of Algorithm 

Modified F-18 Aircraft - Full-Scale Advanced Systems Testbed (FAST). See Pavlock 2011. 

„ . r . ^ „ _ http://hdl.handle.net/2060/20110015358 

Research quality fuel flow sensors installed. 

Airborne Research Test System, 4 th Generation (ARTS IV) 

Nonlinear Dynamic Inversion Contr ol haws. See Miller 2011. http://hdi.handie.net/2n60/2m 1 om s9so 
Autopilots: 


Altitude Hold 
Airspeed Hold 
Wing Leveler 

Algorithm adds biases to: 

symmetric aileron 
trailing-edge flaps (TEF) 
leading-edge flaps (LEF) 


Mode 

Selection 
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Fuel Flow 
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Peak-Seeking 

Algorithm 


Symmetric Aileron, TEF, 
LEF trim position 



Stick/Rudder 
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Output 
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Throttle Command 


Research Fuel Flow Meters 














Technology Maturation (TRL) i 


Technology Transition Map 
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n-D Performance Function 


ID Performance Function 
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Eco Demonstrator? 



Upgrade Existing 
Aircraft 


FAST (F-18 853) 

Multi-Effector Flight Research 
(Prototype in Relevant Environment) 


Single Effector Sim Study on X-48B 

http://hdl.handle.net/2060/20110015999 
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C-17 
State of the Art: 
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Time i 




Flight Research Approach 


Performance Function vs Deflections 

Batch Simulation 

using Simulink Autocode Interface (SAI) 
in f!8sim 



x2: Trailing Edge Flap (deg) 


Performance Function Identification 
Flight Experiment 



Notional Flight Test Point 
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O Baseline aircraft 
O Initial surface biases 
O Algorithm engaged 



Time (approx. 10 minute duration) 


Piloted HIL Simulation 
using ARTS 



Batch Simulation (SAI) 
Using PFI Surface Fit & Noise 
for Tuning and V&V 
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Flight Research Approach 



t Test Point 


O Baseline aircraft 
O Initial surface biases 
O Algorithm engaged 
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Piloted HIL Simulation 
using ARTS 




Algorithm-Engaged Initial 
Flight Experiment 





Fix, Modify, Tune, 
Regression Test 
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Algorithm-Engaged 
Flight Experiment 





Concept of Operations 

Test Flow: 

- Get on condition (25k ft, 240 kcas) 

— Pilot selects mode (configures experiment) 

- Arm ARTS 

- Engage ARTS (autopilots activated) 

- Inserts initial trim biases 

- Turn algorithm on 

- Iterate algorithm repeatedly 

- Turn off algorithm and re-insert initial trim 

- Disengage 



ARTS Engaged 
(Autopilot & Autothrottle) 



Completed Flights as of October 1, 2012 

• Open Loop: Performance Function Identification (PFI) 

- Flight 132, Aug. 7, 2012 

• Autopilot evaluation 

• Matrices of surface deflection combinations (6 completed) 



• Algorithm Engaged: 

— Flight 133, Sept. 19, 2012 

• Ailerons/TEF mode (2D), started from high-drag initial combinations. 

- 3 runs completed 

• Ailerons/TEF/LEF mode (3D), started from high-drag initial combination. 

- 1 run completed 

- Flight 134, Sept. 25, 2012 

• Ailerons/TEF/LEF mode (3D), started from production trim combination. 

- 3 runs completed 

• Ailerons/TEF/LEF mode (3D), started from high-drag initial combination 

- 1 run completed 


Flight 132: PFI Flight Data Examples 




Resume test from H 



_^L 
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18:33:20 


18:36:40 


18:40:00 


18:43:20 
UTC Time 


18:46:40 


18:50:00 


18:53:20 


Delta Fuel Flow 


Flight 132: Estimated Performance Function 

• Recognizable shape 

• Substantial gradient relative to noise 


Estimated minimum fuel flow 



Trailing Edge Flaps (deg) Symmetric Aileron (deg) 


Delta Fuel Flow due to Aileron, TEF, and 
LEF Deflections (for simulation) 


Flight 132: Summary of PFI Flight Results 


Questions Before PFI Flight 

Is the approach feasible? 

• The algorithm detects small changes in fuel flow. Noise 
and disturbances may be too large. 

• PFI experiment will quantify the signal/noise ratio. 

Minimum duration dwell-time interval? 

• Short intervals are desired for faster convergence, better 
use of flight time. 

• Short intervals increase the impact of disturbances. 

• PFI experiment will inform the designers' choice of dwell 
time for the algorithm. 

Can autopilot transients be reduced? 

• Short settling times & minimal overshoots are desired for 
faster convergence, better use of flight time. 

• Autopilot evaluation will include 3 autopilot gain sets. 

What is the shape of the performance function? 

• PFI data will be used to choose initial conditions 

• Surface fit to PFI data will be used in control room to 
verify algorithm is ‘on course’. 

• PFI data will be used in post-flight analysis & technical 
reports. 


Answers from Post-PFI Analysis 

The approach is feasible. 

• Substantial gradients were seen between trim 
configurations despite standard deviations of 
around 50 lbs/hr. 

Dwell time intervals should not be fixed. 

• Lesson learned: Manual advance allows 
flexibility for maneuvering. (Pilot's suggestion.] 

• 30 sec is a good minimum dwell time. 

Autopilot performance is good. 

• Nominal gainset was selected. 

• Good sim prediction of autopilot dynamics. 

• Pilot A: "These autopilots are rock-solid on 
condition." 

Second-order polynomial (paraboloid] fits the PFI 

data well. 

• Six initial conditions selected. 

• Performance function added to sim for algorithm 
tuning. 


Pilot-Selectable Algorithm Parameters 

• Number of effectors (2D/3D mode) 

• Initial trim bias, A through F 

• Number of measurements fed to Kalman filter, M=3, 5, or 7 

• Gain, 7 options from low to high 


Fuel Flow (percent) Fue | Row (percent) 


Flight 133 - First Algorithm-Engaged Flight 



2d, IC:C, M:5, gain:-0.068 




2d, IC:B, M:3, gain:-0.068 



3d, IC:F, M:5, gain:-0.068 
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Flight 133: Algorithm Iterations 
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Note: Fuel flow should be de-trended to account for the airplane getting lighter due to fuel-burn. 

This plot shows the raw results without de-trending. 
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133: Comparison to Estimated PF 




Trajectories versus Estimated Performance Function (Flight Data) 
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Predictions & Results 


Simulation Flight 

Trajectories versus Estimated Performance Function (Simulation) Trajectories versus Estimated Performance Function (Flight Data) 
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- Second Algorithm-Engaged Flight 

Trajectories versus Estimated Performance Function (Flight Data) 
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Flight 134: Detail 

Trajectories versus Estimated Performance Function (Flight Data) 
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Comparison of Completed Flights 
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Near Term: Next Steps (now) 

• Improve fine-tuning performance 

- Tune V and W covariance matrices in KF 

- Longer dwell times, more pre-filtering 

- More measurements (M) fed to Kalman filter 

- Reduce minimum Persistent Excitation (PE) 

• Not a point-design 

- Fly slower, e.g. 200 knots (within the envelope approved for experiment) 

- Fly higher and lower 

• Vary the configuration 

- Empty centerline tank, smokewinders 

- Speedbrake 

- Rudder toe-in 

• Does it work with production sensors? (not research-quality) 

- Throttle position, stock fuel flow meters 


Ideas for Future Work 


• Expanded testing with current platform aircraft (FAST F-18): 

- Stores 

- High speed 

- Coordinate with other upcoming experiments (dual experiment flights) 

• Transition to Transport-Class Airplane 

• UAVs (Ikhana, etc.) - requires more automation. 


Questions and Discussion 


BACKUP SLIDES 


Performance Function Identification (PFI) 

Autopilot, Autothrottle Evaluation 
Matrix of Trim Allocations 



DAG: "Dial-A-Gain" mode selected by 
pilot to select experiment configuration 


• 

DAG 14: 3x3x3 

■ 

DAG 15: 4x4 

A 

DAG 25: 8 large deflections 

* 

8 local by DAG 


Algorithm Flight Test Approach 

• All test points at: 

25,000 ft (+/-2000ft) 12 “^ 

240 kts 


• Initial ICP Evaluation at: 

DAG 17: 2D Nominal Gain 
CAT 17: IC-C, 5 Measurements 

• Plan: 

- If DAG 17, CAT 17 is well behaved 

• Continue to evaluate the other 2D IC's 

• (DAG 17, CAT’s 15-18) 

- If DAG 17, CAT 17 is not well behaved 



oj 

-15 


• Determine if behavior is related to gain or 
number of measurements(M) 


• Select DAG, CAT to attempt to improve 
behavior (different gain and number of 
measurements(M)) 

• Further ICP Evaluations (as time 
allows) 

- 3D Evaluations with same gain and number of 
measurements(M) as successful 2D 

- Evaluations with different gain sets and 
measurements(M) for both 2D and 3D 
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Pilots' Comments 


• I did notice some low freq movement on the first ICP flight. It 
wasn't bad and eventually only noticed when I "looked for it". 
And when I "looked for it" I didn't always find it. 


